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ABSTRACT 

Landscape modeling is the process of overlay and analysis of a large number of multiple 
spatial variables extracted from maps, airphotos, LANDSAT imagery and other remotely sensed 
data. Its applicat-ion to this problem has employed a data base containing 103 spatially re- 
gistered variables. These variables have been assembled to model a site of 292.2 square kilo 
meters at a 1 hectare resolution in the tropical forests of Northern Thailand which are sub- 
jected to the dynamic practice of 
changes in land cover measured with 
terpretation maps were used to model 
trends in 6 forest and 3 agricul- 
in the forest area. Computer com- 
secutive anniversary date LANDSAT 
completed with this data base to de- 
simple image analysis approach for 
this forest cover. LANDSAT image 
within the context of the landscape 
tion of spatial overlays of ancil- 
cation accuracy for the 9 land 
evaluated (slope, aspect, elevation, and distance to drainage) as well as cultural feature 
overlays (distance to roads and trails, permanent housing, and temporary huts). Stepwise dis- 
criminant analysis was employed as the classifier with 3 types of training sets (rectangular, 
aspect, and grid sampled training sets). Two assumptions of the apriori probabilities of oc- 
currence of the 9 land covers within the site were also tested with this classifier, 1.'*., 
equal apriori probability of occurrence and probabilities proportional to the known o.currence 
of each cover type. 

INTRODUCTION TO LANDSCAPE MODELING CONCEPT 
What is Landscape Modeling? 

Geographic data base and remote sensing technologies are about to be combined into a 
whole wnose sum is greater than its parts. This serendipitous meeting of the separate .techno- 
logies will provide a realistic, dynamic infonriation base and associated analysis techniques 
for practical, environmentally oriented land planning. Landscape modeling is the name proposed 
for this fruitful combination, which organizes and overlays data from existing maps, analysis 
of remote sensing imagery, and tabular data into a single computer framework. Tnis assemblage 
provides a multivariate, multi temporal mathematical model which represents the landscape much 
as a topographic map represents a three-dimensional model of the physical land surface. 

Coupled with this composite of data overlays is a collection of computer techniques which allow 
meaningful simulation of the spatial or map-like behavior of this landscape to natural and man- 
induced alteration and control. The current thrust of landscape modeling is the projection and 
display in a map form of the future landscape which would result from the continuation of cur- 
rent land management practices or the lack thereof. Success in this short-run objective has 
enabled the design of the computer techniques needed to predict various scenarios of antici- 
oated landscape alternatives^ The scenarios evaluated could include diverse objectives, such 
as new forest management practices, new zoning patterns for urban land planning, alternate 
sites for a new power plant, or the impact of siting a new dam. Land management should be sub- 
stantially improved if the future spatial Implications of a contemplated action were modeled 
before any coimitment to a fixed cou*"se of action. 

How Does Landscape Modeling Relate to Remote Sensing? 

Computer analysis of remote sensing imagery is symbiotic with the process of landscape 
modeling. It provides the important current and past land cover inputs to the landscape model. 
Detailed interpretation of low- and high-altitude airphotos acquired on various dates initially 
provided the maps used as quantitative measures of land cover change. LANDSAT remote sensing 
imagery, which has been continuously obtained since 1972, currently provides additional infor- 
mation of land cover change. Accurate computer analysis of land cover with LANDSAT imagery of 
various dates provides a more economical, timely, and direct measure of land change for land- 



shifting cultivation. Recent 
4 dates of sequential airphoto in- 
the overall trends and spatial 
tural land cover types occurring 
pari son of the changes in 5 con- 
images for 1972 through 1977 was 
termine the usefulness of this 
monitoring the dynamic changes in 
classification was also tested 
model to determine the contribu- 
lary data to final map classifi- 
covers. Topographic overlays were 
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scape models and their projections of future land cover patterns. 

The accuracy of the computer Interpretation of the remote sensing Imagery Is, In turn, 
substantially Improved by Including a variety of landscape variables, such as slope and dleva* 
tion. Thus, combining the available remote sensing Imagery together with map Information In 
the landscape model provides a basis for substantial Improvements In both activities. Coinci- 
dent, spatially registered overlays of readily available map Information upon the multlspectral 
Imagery of LANDSAT have provided a basis for marked Improvement In the accuracy of their com- 
puter Interpretation to provide current land cf'.cr maps. These Improved, automatically Inter- 
preted maps are of direct use In environmental pUnning. They also “feedback" directly Into 
the landscape model to provide a timely measurement of past and present dynamic tendencies, for 
change In land cover. 

The bulk of this report will be devoted to the LANDSAT Image analysis aspects of landscape 
modeling. A detailed discussion of land cover modeling has been reported elsewhere* How- 
ever, an overview of Its application to this site will provide Insight Into the dynamics of 
tropical forest cover In areas of shifting cultivation. These dynamics must be dealt with In 
any meaningful remote sensing analyses. 

BACKGROUND OF THE FOREST LANDSCAPE MODEL 
Background of the Site Selection . 

The landscape modeling efi'orts described in this report have been underway for approxi- 
mately 6 years. The Initial and more complex application of the approach has been a model of 
the evolution of natural lands to urban lands In the Denver metropolitan area of Colorado*. At 
the Initiation of that study It was ascertained that a parallel application of the approach 
would be made to examine the evolution of forest land to temporary and permanent agricultural 
land. This application was thought, at that time, to provide a cleaner test of the concept, as 
fewer land cover types were Involved* and one overall political control existed for the area 
rather than the more complex spatial political constraints of U.S. urban areas**. Man-made 
changes do dominate the forest area selected but are more closely coupled with the landscape. 
For example, the clearing of trees for shifting cultivation Is known to be a function of lard 
slope, soil type, and other landscape variables. 

At first appraisal It might appear that such a project could be handicapped by a shortage 
of maps, airphotos and LANDSAT Imagery. However, there was no shortage In the quantity of the 
data available, as the particular site has been and Is being closely studied by various Inter- 
national and Thai agencies. Several key limitations did emerge as the study progressed and 
have Impacted upon the approach and results. 

a Inaccessibility prevented careful field control of the land cover maps derived hy airphoto 
Interpretation. 

• No detailed soils map existed for the site although this Is typical of tropical forest areas 
and must be compensated for by any procedures which are to be generalized. 

• Excellent anniversary LANDSAT Imagery exists for every year to date due to an annual cloud- 

free dry season. Unfortunately, this period has been too short to allow the collation of 
meaningful mul tiseasonal LANDSAT Imagery; however, this Is also typical of tropical forested 
areas. * 

Site Description . 

This tropical forest site was selected to represent the northern region of Thailano. It 
IS situated near the Golden Triangle at the corner Intersection of China, Burma, Laos, and 
Thailand (Fig. 1). It represents mountainous area which was orlglnully heavily forested but 
has been rapidly cleared for shifting cultivation of opium, dry-land lice, tea, and other cash 
crops. It Is typical of a large portion of the world's tropical forests where cut-and-burn 
agriculture Is out of balance with regrowth and the mountain slopes and watersheds are being 
rapidly denuded. 

Construction of the Landscape Model . 

The forest area modeled was enclosed In a rectangle of 432 square kilometers (24 km E-H 
and 18 km N-S) which was cellularized with a resolution of one hectare (approximately 2.5 
acres) yielding 43,200 cells. The closed. Irregularly boundeo study site within this rectangle 
was limited to the area connon to several of the special maps made available from the earlier 
analyses of the site by others. It consisted of an area of 292.9 km* or 29,290 square cells of 
lA*_r?59ll/Si5!!.5I!^ anti -watershed bounded by major water courses or drainages on all 

• Nine forest/agricultural land cover types described the Thailand site, while 24 were re- 
quired In the Denver site. 

** The Denver site Is approximately the same size but Is made up of portions of 5 different 
counties which control the land use and thereby the land cover by using 5 different Inde- 
pendent sets of zoning regulations. 
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sides with a mountain ridge In the center. The landscape model of this site consists to date 
of 103 Image and map variables overlaid such that each constitutes a cellular data plane of 1 
ha resolution covering at least the common area of 29,290 ha (Fig. 2). Utmost care has been 
exerted to Insure that each of these data planes registers upon all other? to the nearest 1 ha 
resolution cell*. Input of those data planes derived from existing maps and airphotos was com- 
pleted entirely b> a manual dot sample method**. Area planes, such as topographic elevation 
and geology, are directly sampled cell-by-cell from 1:50,000 scale maps. Additional derived 
ai'ea planes are computed, such as topographic slope and aspect from elevation". Point feature 
planes, such as the location of temporary huts or permanent dwellings, were Interpreted for 4 
different dates from low altitude, black-and-white a1rphotos\ These-planes were then computed 
Into minimum distance area planes*. Linear features, such as drainage and roads and trails, 
were similarly Interpreted from the airphotos for 4 date^'. They were also converted by com- 
putation Into minimum distance area planes*. Airphoto Interpretation maps were prepared for 4 
different dates for 9 forest and agricultural land cover types and overlaid onto the -wCel via 
the area dot sampling procedure". These 4 land cover maps were all Interpreted by a single In- 
dividual In a consistent fashion. 

MODELING CHANGES IN FOREST LAND COVER 
Visual Display of the Changes In Land Cover . 

The 9 types o^ land cover interpreted on each of the 4 dates of airphotos (Figs. 3, 4, 5, 
and 6) Include the 6 forest types of: 

a dry dlpterocarp forest a mixed deciduous forest with ttak a dry evergreen forest 
a dry dlpterocarp forest with pine a hill evergreen forest a teak plantations 
together with 3 agricultural land covers of: 

a shifting cultivation a Irrigated rice paddles a tea plantations. 

Intercomparison of these land cover data planes of the site on a cell-by-rell basis yields 
jraymaps of change which provide qualitative Insight Into the current Imbalance between forest 
depletion by shifting cultivation (Fig. 7) and forest regeneration upon agricultural abandon- 
ment (Fig. 8). A tabulation of the amount of each land cover for each available date Illus- 
trates the rapid recent change In the use of this area (Fig. 9). Clearly the nominal 20 year 
cycle from forest to shifting cultivation and back to forest has been markedly shortened. A 

subsistence level, dry land, permanent agriculture Is setting In with attendant deleterious en- 

vironmental effects, such as the alteration of the basic hydrologic processes. Erosion of the 
soil Is rapidly Increasing with mart^ed Increase In sediment and water yield accompanied by 
shorter duration hydrographs. Thus the cut-and-burn process Is rapidly approaching the Irre- 
versible point where all significant natural forest regeneration will permanently cease. 

Trend Modeling of the Land Cover . 

TKe assumption that future' changes In land cover can be measured In terms of those which 
occurred In the recent pact allows a simple projection of land cover of the area represented by 
a landscape model*. Cross- tabulation of each pair of consecutive land cover data planes In the 
landscape model (e.g. , 1968 with 1972) provides a probability transition matrix. This matrix 
Is applied In a Markov process to the distribution of land cover recorded In the second c.‘ more 

recent land cover data pl.tne to project future trends**. Only a few years of validity may be 

assumed for these projections, as the processes controlling the land cover during the training 
period (I.e., 1968 to 1972) can only be assumed to persist for a few years Into the future 
(I.e., >1972). The particular land cover trends based upon the most recent pair of lano use 
planes In this landscape model predict that Irrigated rice lands will continue to Increase In 
the future as they did during the recent past or training period. However, this cannot be the 
case as the area of land suitable for Irrigation Is being rapidly exhausted even though It was 
not yet In short supply during the training period. This longer term breakdown of the Markov 
trend model points out that It has used only the actual changes In the f'ecent land covers to 
project future trends. However, additional landscape paranieters are available In the landscape 
model to Improve the land cover modeling procedures. 

Spatial Modeling of the Land toyer . 

A much more compel ex mo(ieT has been developed which Is capable of projecting the future 
spatial behavior of the land cover of the site based upon controlling landscape features. Es- 
senTTaTly, It takes all the specific changes In land cover tabulated betv*een the two consecu- 
tive )4*'9_cover_data planes, groups them together Into like changes, and correlates their 

* Technical details on the construction of the data planes whose specific use Is reported 
herein occurs as part of the caption for the figures Illustrating those data planes and 
In the references noted. Illustrations of the data planes (Figs. 3 through 38) have been 
reproduced at a very small size In order to Include the maximum amount of material In the 
limited pages allowed. 

** Land cover trend models have been computed for the following pairs of dates: 1954 with 

1972, 1954 with 1966, 1966 with 1968, and 1968 with 1972*. 
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(*ra<r*-T?nv-» the corryspv^ndlng Und>v.i;'e par.v'vt»'rs In the landscape wdel (e, 9 ., slope, 
distance to ivacls and trails, etc.) usInQ a Multivariate approach. Evacti> the sa^v coiiputa* 
tiona' apiwach Is taKen as in the application of stepwise discvimlnant anal>sis to invioe clas- 
sification except the class to Oe piTdlcted is a specific Chanoe in land cover rather than a 
particular land cover . The variaPles used in this ni»w t>pe of classification are restricted to 
landscape par.i-vti'»*s ano no imaoe param,*ters are included. The needed ceef»ulational procedures 
have Peon worced out .md spatial nvip projections of futui'e land cover have been ccunpleted al- 
though thev ai'd not included here for the sav.e of Previty*. It is contemplated that the cell- 
Py-cell overlay of these spatial "Viul ivesults onto the LANL'NaT images already in the landsvape 
•x'del will prv'vide a Oasis for improved image classification. Apnon kn,H,'edge (i.t., proP.io- 
ilityl of the land cover occurring in each individual image cell based on such spatial p*vje. - 
tions using earlier known land covers and landscape parann^ters should improve the single cell 
classification accuracy. Correspondingly , apriori kni>*ltdde of the overall land cover c.vmposi- 
tion of the site will be subseguently shc«n to improve overall classification accuracv . 

LAMl'bAT OIFFEKENCING FOR F0RE5’ i AN: CClVtR CWAV.i SAPs 

The simplest LAN?SAT image analysis undertaken was the cell-by-cell subtraction of sach 
ICS band of one image frvm those of anotner ann’ve»'sar> imaee, Th« radiance differences which 
result rang# i about jero and the i magnitude of these differences correlate with various 
changes in land cover between the two images. For example, clearing the forest or brush cover 
to expose the soil and oead g*vund litter for shifting cultivation will increase me radiance 
in iCS-b and decrease it in MSS-’, causing image differen.es greater than rero and less than 
:ero, respectively*. Annivei'sarv date image di '♦ei'encing and ratioing has rvt yet been widely 
tested due to the difficulties and costs associated with t^e rectification of several kAhi^SAT 

scenes to a con-xm cellular base. The supplying of already rectified imagery from LMrSAT J by 

NASA tiSFC and the [itOS Data Center shc'uld inc*'ease the interest in these proce.iu*-es and their 
testing. Ground truth collection for the identificatuxi of the changes detected is also cx'>-e 
complicated as it must be perfoirx'd twice. 

Application of LANDSAT image differencing to mar chan.;es in t»\>p1cal forest cover was 
testec on this site since an understanding of its dynamics was available. The site is also 
small so that independent gev'metric r^.-ti f <cat ipr w.is .feasible for e.ich lAhi’^AT anniversary 
Image of 19’?, 19’’A, 19’S, 19’o, and 19"'’. The 4 o.inds representing each o* thexe images 
were first transfc’mrd to remove the major sensor distortiv'ns and resancled into the 1 ha 
square Cfll data planes using a nearest neighbor algorithm*. This new cell is approximately C 
tii'v's the area of the original LAHPSAT picture elei*v’nt and thus about bOt of the LANCSAT ele- 
rmnts were transferred to this new format (Figs. 11 and 1^1. The topv>graphic and otner land- 
scape maps introduced into the landscape model wei-e accurate only to the 1 ha cell $ise and the 

LAN.'^AT image>> was undersanpled to match these data planes. This unde*-sampling procedure has 
not appeaised to introduce any noticeable effects on any of the LAN.ISAT image processing activi- 
ties except the incorporation of some single cell noise into tn# difference gravewps whprt the 
i independent rectifications do not sample precisely the same nea'^est neighbor. The effect of 
most of this noise was ror\'ved by the application oV a Km band pass spatial filter to each 
difference image befote its display as a graymap*. 

little experience with image differencing can be found in the literature in connection 
with tropical forests and shifting cultivaticm. Thus the initial question addressed was what 
time base shciuld be used between images le.g., 1 year, : years, etc.'. Subsequently, all com- 
binations were computed and 1 year intervals found to produce the nc*st cons'Stent results. The 
image difference data planes have a narrow histogram of Gaussian appearance in all cases if th» 
$pl.ir angle differences between the 2 dates is no more than a few degrees. Location of the 
data ranges to either side of the mean which contain the land cover changes sought is not 
straightfoiward unless accurate ground cv'ntrol is available for both image dates. The histo- 
gram of the whole image difference data r’ane for a given *»ss oang reverses from positive to 

negative if the order of subtraction of the two images is reversed. This Is the sa:^ as chang- 

Ing the arithmetic sign of the solar angle diffe»ence. hv'wever, this arithmetic sign can also 
be changed by a small change of 1 or ? weeks in the exact calendar date of one of the images 

selected to represent the “anniversary*' dates. This also causes a positive to negative reverse 

of the histogram of this data plane. Considerable additional experience with this simple pro- 
cedure is still needed to obtain a moii# generali;ed understanding of its application to ttimi- 
toring forest cover type and change. For examt'le, inc»-easing the solar angle difference be 
tween the two images used enhances the contrast of one Forest Avpe relative to another m the 
difference image. This might be due to the inherent differences in the fee stand tixi»*pho’ogv 


• The mean values of these difference data planes can fluctuate several units eithe*- way ♦npm 
jerp due to overall image differences between the two anniversary dates such as those in- 
duced by different atnvspheric abso*Ttic*n in one image versus the other. 
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of the 2 forest types which cause different subtle shadiw* effects in the images selected, thus 
alK^ing discrimination of these types in the difference graymap. 


Consistent detection of the area of shifting cultivation was observed for each of the 
graymaps produced to display the 1 year anniversary difference images of mss- 7 (figs. 13, 14, 
15, and 16). The first two or earlier date difference gravm.ips enumasice the general pattern 
of Shifting cultivation (Figs. 13 and 14). The remaining two wre recent graymaps (Mgs. 15 
and 16) ens’hasije the active areas of shifting cultivation at the hignest altitudes at the 
center of the site which have been cleared and planted to opium popple-, during the intervals 
between images. The relative detection of old versus active shifting cultivation in this se- 
quence of 4 consecutive 1 year interval diffei-ence graymaps is not clearly understood although 
it IS due in part to the nanwwness of the histograms noted earlier which markedly restricts 
the exact selection of the difference data ranges to be displayed in the graymap. 

Consistent patterns of the imvvement of permanent agriculture into the foi'ested portion of 
the site can be detected annually in the MSS-5 image difference graymaps (Figs. 17, 18, 19, and 
20). Unfortunately, these changes can be confused with radiance differences in the areas of 
agriculture present on both dates where changes in crop type or planting date occur from year 
to year*. However, careful analysis of the original airphoto land cover maps shivws a steady 
progression of permanent agriculture into the foivst from the NE and SE corners*. This steady 
encroachment of irrigated and dry land r*ce, together with orchards, can readily be detected 
on the difference graymaps. 


Image differencing also tends to compensate in part for topograpnic induced radiance vari- 
ation if the solar angle difference between the images is not too large. Mi>re meaningful com- 
pensation for these topographic effects would result if the cell-to-cell image diffei'ence were 
weighted by the radiance value of each cell for one of the input images. This “automatic gain 
control" would adjust the image differences tor similar land Cviver changes in the shadv^ws or 
sunlit areas to the same numeric value. 


LANOSAT forest land cover classification with ANCILLARi DATA 
Classification with Topographic Overlays . 

LANOsAT classification in areas of significant relief is limited in accuracy by that por- 
tion of the radiance variation created by the terrain which is inde^'ende'lt of the surface cover 
type. Deep shadows represent the limiting condition but surface areas which are not in shadows 
out merely at different topographic slopes and aspects can exhibit significant topographic in- 
duced radiance variation within a specific land cover type. A large portion of the economical- 
ly important forested areas of the world occur in areas of mOiierate to severe relief and these 
are also the areas of least knowledge of the forest type, condition, and product! vity . The 
forest site selected for this analysis is typical of the rugged tropical forest remaining in 
Thailand and adjacent Asian nations. It is the very ruggedness of the terrain which has pre- 
served a portion of these natural forests to the present time. A quick, comparison of the rec- 
tified 1973 LANDSAT images (Figs. 11 and 12) with the 1972 land cover clearly illustrates that 
the radiance distributions Imaged by LANDSAT are dominated by topographic effects and shiw» 
little resemblance to the land cover (Fig. 6)**. The procedures developed for *he application 
of LANOSAT classification methods to such areas must cope with these terrain induced varia- 
tions. Three different attacks on this problem have been initiated as follows: 

• Selection of training sets to represent each cover type as more than one class according 
to its slope and aspect range. 

e The overlay of topographic features as data planes onto the imagery and their incorporation 
as input variables in the classification. 

e A deterministic adjustment or preprocessing of the LANDSAT MSS bands to remove the major 
portion of the topographic induced variation. 

The results of the first two approaches are reviewed while the third has not yet been success- 
fully completed although a graymap of the site ve>^ similar in appearance to that of MSS-7 
(Fig. 12) has been simulated . This conf'utation of the radiance variation in MSS-7 uses the 
topographic slope and aspect data planes and other data to compute a new data plane consisting 
of the cell-by-cell near- instantaneous solar insolation at the exact time of the LANDSAT im- 
age*. 

The topographic overlays employed in the two empirical approaches to terrain compensation 
were derived from the 1:50,000 topographic map of the site. The elevation was sampled for each 
x'ells yielding the initial data plane (Fig. 21). This cell by cell deter- 

* This is not necessarily an undesirable conditiew, but one we do not understand how to sort 
out and use to advantage. 

•• Note that the December airphotos used in preparing the 197? land cover map are separated by 
about 1 month from the 27 January 1973 LANDSAT image and thus represent the same year and 
dry season. 
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niln^tion Mas vpr> tedious but yielded the most accurate elevation data plane of any technique 
available*. The slope (Fig. 22) and aspect (Fig. 23) data planes Me**e obtained by fitting a 
faceted surface to the elevation data plane and computing the slope and aspect of this surface 
at each of the 29.290 cells. The drainage iwp Mas Interpreted from the 1072 airphotos and con- 
verted to a binary data plane. Each remaining cell not occupied by a drainage element Mas as- 
signed a value equal to the nunt'er of cell units (1.#., distance) to the nearest cell occupied 
by the drainage In a N-S and or E-W excursion (Fig. 24). These 4 topographic data planes Mere 
overlaid In a cell-by-cell fashion onto the 4 *tSS bands and 6 ratios of HSS bands In the land- 
sc.ipe model. 

SteiMvise discriminant analysis Mas employed as the classifier In these analyses as It Is 
reasonably forgiving of the non-GauSSlan nature of the data to be classified. Extensive train- 
ing set classification Mith a Mide variety of Invige and landscafH* variables and types of train- 
ing sets Mas coniileted. This provided the basis for the selection of the training sets and 
variables to be tested for their final verification or map classification accuracy. These fi- 
nal verification results m111 be summarised here by the Illustration of the actual classifica- 
tion maps and their verified accuracies. At this point It Is appropriate to note that, based 
upon training set accuracies. It Mas determined that the 6 ratios of the 4 MSS bands contrib- 
uted little to the classification accuracy and they Mere omitted from further analyses. The 
Impact of three types of training sets tested m111 be presented and they are: 
e RECTANGULAR TRAINING SETS consisting of several rectangular groups of cells selected to 
represent each of the 8 cover types** as Identified on the 1972 land cover data plane. 

• ASPECT TRAINING SETS Mhere each of the 8 cover types Identified on the 1972 land cover 
data plane, except Irrigated rice, Mere subdivided Into 2 new classes and represenred by 
rectangular groups of cells consisting of the particular cover type on sun-facing and 
opposite-facing topographic aspects. This yielded 15 classes Mhich Mere regrouped back 
Into the 8 original classes after classification bu before display or verification. 

e GRID SAMPLED TRAINING SETS where every third row and column of cells In the 1972 land 
cover data plane Mere extracted and regrouped Into 8 training sets representing the land 
covers. 

Two assumptions for the apriorl probabilities of occurrence of each land cover were tested as 
follows : 

e EQUAL APRIORl PROBABILITT of occurrence of each cover type In the area to be mapped. 

# proportional APRIORI PROBABILITY of occurrence of earn cover type In the area mapped, 
where the probabilities were taken as directly proportional to the populations of each 
of the cover types In the grid sampled training sets. 

The 3 training approaches, together with the 2 types of Initial assumptions of apriorl 
probability of the final map composition, yielded 8 classification maps with MSS bands 5 and 7 
(Fig. 25), the 4 MSS bands (Fig. 26), and the 4 MSS bands with topographic variables (Figs. 27 
through 32). Visual Insight Into the success of these approaches to mapping shifting cultiva- 
tion can be directly observed from the classification maps. Overall performance of the clas- 
sifiers tested Is best evaluated by the cell-by-cell verification of all 29,290 cells In the 
1972 land cover data plane to be subsequently explained and presented. 

Classification with Cultural Feature Overlays . 

The distribution of three types of cultural features was mapped from each set of airphotos 
and converted Into a data plane, specifically, the location of permanent housing, temporary 
huts, and roads and trails (Fig. 33). Each of these 3 types of features was computed Into min- 
imum distance planes In a fashion similar to that of the drainage plane (Figs. 34, 35, and 36). 
These new data planes were overlaid upon the 4 MSS bands and 4 topographic variables to deter- 
mine their additional contribution to the classification results***. Since there Is a close 
correlation between the location of these cultural features and the surrounding land cover it 
seemed reasonable that they would further improve the classification maps (Fig. 37). 

Verification of the Accuracy of the Classification Maps . 

various classifiers tested was directly evaluated since the 

* Practical future applications might employ something like a cell-by-cell stereoscopic sat- 
ellite capable of providing the synchronized elevation values needed for terrain compen- 
sation during land cover classification of multlspectral Imagery. 

** The ninth cover type, representing teak plantations, occupied only 40 ha of the site and 
was omitted for these classifications. 

*•* This may appear to employ circuitous logic; In other words, that such cultural features are 
part of the desired solution and not part of the Input data. However, we have been careful 
to use the words "land cover" whereas these cultural features represent "land use". They 
might well be quickly mapped as point and line features from airphotos or satellite side- 
looking radar and aci^d to the LANDSAT classification procedure to assist It to "paint" In 
the detailed areas of land cover. 
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".inNwiM-** Is ahTicl) available for each of the ;N,2'>0 cells classified In the form of the 197» 
land cover data plane. A cell-by-cell conisiarlson of the particular classification map with the 
land cover map yielded a simple verification displav where a cell can qulcl^ly be Identified as 
right (black) or wrong (white) If Its classl flcatlo»< checks with the kn^^wn land cover (Fig, 5^). 
Thoc verification displays are useful to quickly ch’temilne the overall results of the classi- 
fication and determine If any spatial Mas exists In the classifier. The percentage of each 
cover type correctly classified, together with the overall accuracy , directly compares the 
v.irious approaches (Table 1). 

It must be noted that the percentage of con-ectly mapped points in such a verification 
cannot be expected to closely approach ioot. The airphoto interpretation maps of land cover 
have theli v>wn inherent accuracy which Is not 1001. This initially Impacts on the training set 
Identification with these land cover maps. The Impact Is even greater during the verification 
of the LANDSAT classification results with the "known’' airphoto derived land cover maps. The 
verification accuracy can only asymptotically appn^ach the accuracy of the airphoto land cover 
map with the Inclusion of each additional variable or set of variables. For example, assunv 
that the land cover map derived from the airphotos Is 70* correct; that Is, If 10 random cells 
wo'e checked on the giwund 7 of the 10 would have been assigned by the photo interpretation to 
the correct cover type*. The classification map may also be '0*. correct relative to ground 
condition*. Thus, verification of the lANOSAT classification map against the airphoto map will 
pivduce a 70S X to*: or 49'» verification accuracy. Whei'e training sets have been identified 
from the samt' airphoto land cover map a tendency will exist to bias the classification solution 
tvwijrd the map results and away from the real ground conditions somewhat Increasing this 49* 
figure. 

The verified map classification achieved the best overall accuracy of 55. 4t with the grid 
sampled training sets** and proportional probabilities. Examination of the results of this ap- 
proach for each combination of va'^lables coim^uted Illustrates that It does asymptotically ap- 
proach the limit set by the accuracy of the ground "truth" In the form of the 1^7r land cover 
map (Table 1). Independent tests of grid sampled training set against other training set ap- 
pixiaches for the Penver meti'opolltan landscape nv'del also pnvduced similar superior results 
even though a lai'ger nunt>er of classes were sought*. This grid sample approach does not pro- 
duce as Mgh a training set accuracy, as It has been deliberately designed to Include all the 
spectral variability In the area to be mapped. Rectangular or other area type training sets 
;w»rform best when training set accuracy Is the criterion In that they have been deliberately 
selected to minimize spectral variability. Subsequent application of rectangular or areal 
training sets to the "real" world, I.e., the area to be mapped with all Us spectral diversity, 
has demonstrated that high training set accuracies do not extend into higher map verification 
accuracies*. 

Training set accuracies achieved with grid sallied training sets are essentially equal to 
the final man verification accuracies since the grid training sets are selected to be repre- 
sentative of tho total area to be mapped. This provides one additional Important advantage not 
possessed by rectangular or areal training sets. Ar\v optimization done upon grid sampled 
training sets will also optimize the map classification operations, for example, the selection 
of an optimal subset of spectral bands for computational econoi^ from among a suite of multl- 
spectraUmul t1 season Images. Optimization of these features based upon rectangular or areal 
training sets does not necessarily extend to the final map classification employing these 
training sets. 

Obtaining grid or random single cell training sets may appear loo tedious to be practical, 
but this Is not the case. First, the sample need not be nearly so large a percentage of the 
area to be mapped as was the case In this research (I.e., 11. IS). Secord, the approach In- 
creases the accuracy of the photo Interpretation and »*educes the human tedium and monetary 
cost of acquiring the training data. All the photo Interpreter has to do Is accurately Iden- 
tify the land cover at specific preselected points on Individual airphotos or stereo pain. 
Thei'e 1$ no need for the error prone spatial generalization which Is required In making maps or 
to a lesser extent In defining rectangular or Irregularly shaped areal training sets. 

* A Treasonable assumption In this case, as seen from the large amount of spatial generaliza- 
tion which exists in the 4 airphoto land cover maps. 

*• The grid sampled training sets used here are egulvalent to random samples taken frxrm the 
ii.TK' arcs If th* sample sizes ai-e sufficiently large. Using a grid Is simply a convenient 
way to extract tht samirl** from oellularlred data. 
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Figure 1. Location of t>c Sticy Site*. 

The irregular study site is 292 <quare kilometers in Northern Thailand near the area of the in- 
tersectiOii uf tne China, Burma, Laos and Thailand borders called the Golden Triangle. It is 
representati ve of the genera' area which is mountainous and was originally heavily forested but 
has been rapidly cleared for shifting cultivation of opium, dry>land rice and other cash crops. 
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Figure J. Conceptual iiatky^ cf t>€ Um'scape -I'icl Constructed for niE Stud> -\rea‘. 

Listed are the 103 variables which have been overla'd as Individual data planes In t*ie fonn of 
:9,:90 sgua»-e cells of 1 ha resolution. Each of these data planes contains 43. .'DO cells with- 
in the Inscribing rectangular area but the 1nro»*mat1on content Is limited to the iri'egularly 
bounded area corrvn to several of the special maps m.ide available bv other studies of f'e area 
This Irregular site consists of an anti-watershed bounded bv watercourses or dra”'av:es 

on all sides with a mountain ridge in the center. 
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Figure 3. 1954 Land Cover Type Map*. 

tlaok ^akifting sultioation and t*a 

nediun gvay-m^ixvigatad rice 
light gvccy—m^hill avsrgreeri foraat 
vhita ■ ^ 4 other foreat typea 
A land cover map was interpreted from 19 medium altitude 
1:60,000 black and white air photos of Jan. 1954 and com- 
piled at 1:50,000. The map was, in turn, digitized to 
record the spatial distribution of each of the 8 land cover 
types in the 1 ha cellular data plane format. 



Figure 4. 1966 Land Cover Type Map*. 

blaok ^shifting aultivatio>t and tea 

medium gt>ay~^ irrigated rise 
light gray—^Hll evergreen foreat 
white ^ 4 other foreat typea 

A land cover map was interpreted from 76 low altitude 
1:20,000 black and white air photos of Jan. 1966 and com- 
piled at 1:50,000. The map was, in turn, digitized to 
record the spatial distribution of each of the 8 land cover 
types in the 1 ha cellular data plane format. 



Figure 5. 1968 Land Cover Type Map*. 

black ■ ‘^ahifting cultivation and tea 
medium grccj -^irrigated rice 
light gray—^hill evergreen foreat 
whi te ■■ -^4 other foreat typea 

A land cover map was interpreted from 70 low altitude 
1:25,000 black and white air photos of Jan. 1968 and com- 
piled at 1:50,000. The map was, in turn, digitized to 
record the spatial distribution of each of the 8 land cover 
types in the 1 ha cellular data plane format. 



Figure 6. 1972 Land Cover Type fiAP*. 

blaaJi ^ehifting cultivation and tea 

medium gray -^irrigated rice 
light gray—^hill evergreen foreat 
white' ^4 other foreat typea 

A land cover map was interpreted from 75 low altitude 
1:20,000 black and white air photos of Dec. 1972 and com- 
piled at 1:50,000. The map was, in turn, digitized to 
record the spatial distribution of each of the 8 land cover 
types in the 1 ha cellular data plane format. 
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Figure 7. 1968 to 1972 f^si [fepirrioN 1V‘. 

A'ni'^rTsJ tc' ahiftin^ juttiVsXtijn 
i::rk ocnv^rt^d tc t4a pl<sntati 

»“/wi’<F#t .vni'jrt^i tw’ :rr:’y’.:r^i jn^ paiJiVa 
Obtained by a ceiwby-cell comparison of the 1972 land cover 
(FI9. 6) with that of 1968 (Fig. 5). 
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Figure 8. 1968 to 1972 Forest Regeneration IV*. 

>: j«( — jst.xndc*'^d la*tJ r^t^twing tc fcx^'ctc 

d.zrk ♦■nwu’ tea plar.tatis:".s 

''i'rKcat4.i rt\v pcJdi** 

Obtained by a cell -by-cell comparison of the 1972 land cove»* 
(Fig. 6) with that of 1968 (Fig. 5). 
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Figure 9. Historical Forest (bvtR 
DEPirrioN AND T>€ (Concomitant Increase 
IN T>€ /V?ea of Shifting Cultivation. 

The area of each land cover type was 
obtained by a simple computer tabula- 
tion of the 1 ha cellular land cover 
data planes for the dates available 
(Figs. 3, 4, 5, and 6) . 
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Figure 10. Predictions of Future Trends 
IN THE -VujNT OF Ea*:h Forest Ccver Type 
AND G>feting Agricultural Unt Use*. 

Simulated usinq a Varkov process and 
transition matriii provided by the C9.C90 
cell-to-cell comparisons of the 1968 to 
1972 land cove** data planes vf^i9s* 5 ana 
6). Similar simulations have been com- 
puted using the transition matrices pro- 
vided by each of the consecutive pairs 
of dates of 1954 to 1966 and 1966 to 
1968 as well as the longer overall time 
interval of 1954 to 1972. These models 
illustrate the tendencies for land cover 
IRRIGATED RICE changes in the study area due to evolv- 
\ ing forest manager^nt policies. 
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Figure 11. imAJ *SS-5 Grav>w> . 

-w-so. :.Y A 

This LANOSAT red (0.6 to 0.7 um) image was ta^en Jan. 19'3 
and has been rectified and resamoled using a nearest neighbor 
algorithm to match the 1 ha cellular jata plane format. 



Figure 12. LAISAT ‘^-2 Gf^A>^wF', 

^ 'vro'." :.v fh.'v 

This LANOSAT photo infrared (0.8 to 1.0 v.ml image was taKen 
Jan. 1973 and has been i*ectified and resa;npled using a 
nearest neighbor algorithm to rnatci the 1 ha cellular datJ 
plane format. 
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Figure 13. 1974-1973 LANDSAT ftS-7 Difference (jRAVMAP*, 
(SMLUisr ss‘]at:vs ra::a,\cs :k suck' 

The Inwqes were t.iKen 22 Jan. 1974 .ind ?7 J.in. 1973 yielding 
a lolar angle difference of -1*. Each Image was Independ- 
ently rectified and resampled with a nearest neighbor algo^ 
rithm tn match the 1 ha cellular data plane format. The 
1973 Image was subtracted from the 1974 Image on a cell-by- 
cell basis. A low frequency spatial band pass filter was 
then applied to remove most of the single cell noise. 
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Figure 14. 1973-1974 LANDSAT MSS-7 Difference fjRAYMAP*. 

mptavcs ^tffffsscz'^ rx buck' 

The Images were taken 13 Feb. 1975 and 22 Jan. 1974 yielding 
a solar angle difference of ♦!*. The difference Image was 
computed just as In Figure 13 above. 
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Fl -AJRE 15. 1976-1975 lANTCy\T i^-7 DlFFEkKNCE liRAVMAr*, 


*»**♦#. •• » • 


UIMEFT r'Sinvr HAC.-^.nr riFFrsFSCyF :x KACKi 

The Images were taken Feb. l'i7o .md 13 Feb. 1975 yielding 
solar angle diffe»"ence of ♦.r'. The d1ffe»‘ence image was 
f o.Mnputed Just as in Figure 13 above. 
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Figure 17, 1974-197^ L/WSAT MSS-5 Difference Graymap', 

{ZAFCKsr pofrrm MciA^tcr crrnFSfiCKs nt hack) 

Tht Inwigcs M«rt taktn 22 Jan. 1974 and 27 Jan. 1973 y1«1d1nq 
a solar an^U dlff*r*nco of -1*. Each Imaq# was Indaptnd- 
•ntly r«ctlfl#d and r*sampled with a n*ar*st neighbor alqo- 
rithm to match the 1 ha ctllular data plana format. Tht 1973 
Imaw was subtracted from tht 1974 Ima^t on a ctll-by-cfll 
basis. A low frequency spatial band pass filter was then 
applied to remove most of tht single cell noise, 


Figure 18, 1975-1974 LAHKM MSS-5 Difference Graymat'. 

iSMALLSST sSuATnr rav:as\'s rrFFKSKWFs :.v 
The Images were taken 13 Feb. 1975 and ?2 Jan. 1974 yielding 
a solar angle difference of The difference Image was 

computed Just as In Figure 17 above. 



Fuuhe B, 197G-1975 

(SMAir.SST ST'^TTXr 


bVflSAT MSS-5 Difference ijRAMw'. 

SAriAXi'F rTFFFFFT'FF :.Y 


The Images were taken 2b Feb. 197b and 13 Feb. 1975 yielding 
a solar angle difference of >3*. The difference Image w.ss 
computed Just as In Figure 17 above. 



FiciURE 20, 1977-1976 LANDSAl 1^-5 Difference liRAYmp', 

■ ^mLLSFT SFOAT 'tT SAC'lAXrF rTFFFFKS,'SF :.V 5UCA ' 

The 1ma,jes were ta*,en 2 Feb. 1977 .snd 26 Feb. 1976 yielding 
a solar angle difference of -7*. The difference Image was 
computed just as In Figure 17 above. 
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Figure 21. lofxx-rfWHic Flcvatu^ 

iSTJHSir KllVATTOSS IS SIACK) 

The elevation of each Individual 1 ha cell wa^ eetlmated to 
el thin 10 meters trow a 1. SO, 000 iCale topographic map for 
each of the cells using a uot grid overlay. The con- 

tour Interval between gray steps of the display Is X)4 m 
(1 ,000 ft). 



Figure 22. Tcvograpmic Scope V*. 

^STSSFSS‘^ s<LCPSii IS HA.'K) 

The slope of each Individual I ha cell was contputed by fit- 
ting a regi'esslon least- square plane to the topographic 
elevations of a 3 by 3 array of ceMs and determlnlna the 
slope of that surface to the nearest 0.1% and assigning It to 
the center cell of the array. Repeating the process over the 
whole elevation plane provides this slope data plane. 



Figure 23. Topuj^aphic Asf*ECT V‘. 

The aspect of each Individual 1 ha cell was computed by fit- 
ting a regression least-square plane to the topographi. 
elevations of a 3 by 3 array of cells and determining ti • 
aspect (1,e., cimn'ass direction of the perpendicular to thv 
plane) of that surface to the nearest \" and assigning It to 
the center >*^11 of the array. Repeating the process over the 
whole elevatlT^TNqJane provides this aspect data plane. 
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A drainage uvap was Interproted from 75 of 
blat^ and white air photos of Dec. 14 ?a. 
turn, uigitlied to inecord the presence or 
age in the ' ha cellular data plane foimiat. ^ 
tance algorithm was aoolled to this binary data 
compute the minimum distance fixim each erx>ty cel 
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pied by the drainage to the nearest cell through which the 
drainage paisds. 
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FiaiRE 25. 2 -Variable Class ificatutj Mv' 


j'hi t* 4 '' . r/'iwr tj,r 

The inwije was taken ?7 Jan. 1973. The 8 cover types we?*e 
classified usiny training sets derived from a 3 by 3 gnd 
sami'le of the 197J land cover data plane. The 2 variables 
analyzed were MSS-5 and MSS-7. The cl assi f Icatlon employed 
ste^>w1se discriminant analysis with apriorl probabilities 
dltoctly proportional to the known occurrence of each land 
cover type In the site (I.e., pioportlonal to the popula- 
tions of these training sets). 


Fu'^e 26, ^-Variable Classificatkti ^Uf 
. samplft r.w.^Y^.v • ff?s \:;r 
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i. un j7»* ^ r er»r .j.x reZ rt 

v>t;e^ .'r*’vr /V-rvet tXtV'tsa 

This classification was Identical to that noted In Figure 25 
above except that It operated upon the 4 variables of MSS-4, 


Fk'AJRE 27, 8 -Variable Classification ^Iap 


i*:\v 

li^ht /''rvaf 

•jht f- • t' .'thtir .'Vivat 

This classi ficatlon was Identical to Chat noted In Figure 25 
above except that It operated upon the 8 variables of MSS-4, 
5. 6, and 7, slope, aspect, elevation, and distance to 
drainage. 


Future 28, S-Varl\bl£ Classi fi cat ktj 'W’ 
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This classification was Identical to that noted In Figure 27 
above except that It enH'loved equal apriorl probabi 1 1 ties of 
occurrence of each of the .8 cover types. 
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Figure 33. Cultural Features IV’. 

: SHOWS SITE BOUNDARY, ROADS AND TRAILS, TST.TJRARY HUTS, 

AND PER^NENT HOUSING FEATURES) 

A cultural features map was Interpreted from the 75 low alti- 
tude black and white air photos of Dec. 1972. The map was, 
in turn, digitized to record the presence or absence of each 
of the cultural features in the 1 ha cellular data plane 
format. TEMPORARY HUTS represents individual brush 5hel ters 
used for temporary or seasonal occupancy. PERMANENT HOUSING 
represents schools, field station, tea processing and storage 
areas, permanent dwellings, etc. ROmOS AND TRAILS range any- 
where from a heavily used footpath to improved dirt roads. 


Figure 34. V of Miniimum Distance to Permanent Housing". 

(SHORTEST DISTANCES IN BLACK) 

A binary data plane was extracted from the cultural features 
data plane (Fig. 33) recording the presence or absence of 
areas of permanent housing in each 1 ha cell. The minimum 
distance algorithm was applied to compute the minimum dis- 
tance from each empty cell not occupied by features of per- 
manent housing to the nearest occupied cell. Diamond shaped 
patterns result from restriction of the minimum distance com- 
putation to combinations of N-S and E-W distances for reduced 
computation costs. 





Figure 35. V of Minimum Distance to Temporary Huts". 

(SHORTEST DISTANCES IN BLACK) 

A binary data plane was extracted from the cultural features 
data plane (Fig. 33) recording the presence or absence of 
temporary huts in each 1 ha cell. The minimum distance 
algorithm was applied to compute the minimum distance from 
each empty cell not occupied by a temporary hut to the near- 
est occupied cell. Diamond shaped patterns result from 
restriction of the minimum distance computation to combina- 
tions of N-S and E-W distances for reduced computation costs. 



Figure 36. V of Minimum Distance to f^JAOS and Trails". 

(SHORTEST DISTANCES TN BUCK) 

A binary data plane was extracted from the cultural features 
data plane (Fig. 33) recording the presence or absence of a 
road or trail in each 1 ha cell. The minimum distance algo- 
rithm was applied to compute the minimum distance f’^om each 
empty cell not occupied by a road or trail to the nearest 
occupied cell. Diamond shaped patterns result from restric- 
tion of the minimum distance computation to combinations of 
N-S and E-W distances for reduced computation- costs . 


13 




Figure 37. 11-Variable Classification ^Vvp^ 

(JRW SAtfFLEV TRAISZltii SFTS AXZ PFCPCRTIOML 
APRTCRI PRCPABrLrrZF:?) 

bljtyk ouZf-t jr.d 

-wiu4r r!*..v 

It^ht ^hi’l 

tw ^ 4 ^'thgv 

The 8 land cover types were classified using training sets 
derived from a 3 by 3 grid sample of the 1972 land cover data 
plane. The 11 variables analysed »<ere MSS-4, 5, 6, and 7, 
slope, aspect, elevation, distance to drainage, roads and 
trails, permanent housing, and temporary huts. The classifi- 
cation employed stepwise discriminant analysis with apriori 
probabilities directly proportional to the known occurrence 
of each land cover type in the site (i.e., proportional to 
the populations of these training sets). 



Figure 38, 11-Variable Verification ^Up^ 

(GRIP SA.^!PLFP TRAIMRG SF7S AXC FRGPGRTIGXAL APRICRI 

PR03ABTLI7TES - HLAFK PEXOTSS PFIXTF GrRRF'TlY 'LiGGlFIFP' 
This display shows spatially those points which were correct- 
ly identified in the 11 variable classification map of Figure 
37. Each and every 1 ha image cell has been classified into 
1 of the 8 land cover classes and then checked against the 
corresponding 1 ha cell in the 1972 land cover data plane. 
When a 1 to 1 match occurred a black cell was displayed. 


ORir.TNAI. PA(;K is 
OK Ol’M n*>‘ 


Table 1, Ofarison of Verification Resilts for t>€ L/VIDSAT Forest Cov^r Cuvssifi cations. 

The numbers represent the percentage of the total number of the 29,290, 1 ha cells which were 
correctly classified into their respective cover type when checked against the 1972 land cover 
map. The second set of percentages in italics represent only those cells correctly classified 
into shifting cultivation. These italicized values for shifting cultivation fluctuate widely 
as the various classification approaches were optimized to achieve the best overall accuracy 
for all 8 land covers and not optimized for this specific cover type. 

2 variables used ■ MSS-S and 7, 

4 variables used • those above plus MSS-4 and 6, 

8 variables used ■ those above plus topographic slope, aspect, elevation, and distance to 
drainage , 

11 variables used ■ those above plus distance to roads and trails, permanent housing, and 

ten^orary housing. 
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